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Lossless filtering of a single coherent (Schmidt) mode from spatially multimode radiation is a 
problem crucial for optics in general and for quantum optics in particular. It becomes especially 
important in the case of nonclassical light that is fragile to optical losses. An example is bright 
squeezed vacuum generated via high-gain parametric down conversion or four-wave mixing. Its 
highly multiphoton and multimode structure offers a huge increase in the information capacity 
provided that each mode can be addressed separately. However, the nonclassical signature of bright 
squeezed vacuum, photon-number correlations, are highly susceptible to losses. Here we demonstrate 
lossless filtering of a single spatial Schmidt mode by projecting the spatial spectrum of bright 
squeezed vacuum on the eigenmode of a single-mode hber. Moreover, we show that the first Schmidt 
mode can be captured by simply maximizing the fiber-coupled intensity. Importantly, the projection 
operation does not affect the targeted mode and leaves it usable for further applications. 

PACS numbers: 42.65.Lm, 42.50.Dv, 42.50.Ar, 42.65.Yj 


I. INTRODUCTION 

Sources with a perfectly single-mode spatial spectrum 
are desirable but rare - an example is a laser with the 
beam quality factor = I. Most sources contain mul¬ 
tiple spatial modes, which leads to the need for filtering 
methods. Preferably, these methods should be lossless, 
i.e., they should maintain all energy contained in the fil¬ 
tered mode. This is important for laser sources but be¬ 
comes absolutely crucial for certain types of nonclassical 
light, because of the destructive role of losses. 

Although formally one can choose free-space modes in 
many different ways, the most common example being 
plane-wave modes, spatial coherence dictates a special 
choice of eigenmodes for each type of radiation. For in¬ 
stance, if a mode has to contain all spatially coherent 
radiation, it has to be chosen according to the Mercer 
expansion [1] of the first-order Glauber’s correlation func¬ 
tion. The Mercer expansion provides the so-called coher¬ 
ent modes because the radiation within each of them is 
coherent. 

Very similar to Mercer expansion is the Schmidt de¬ 
composition. While the Mercer expansion describes co¬ 
herence of partially coherent light, the Schmidt decom¬ 
position also accounts for photon-number correlations. 
Within a certain Schmidt mode the radiation is co¬ 
herent and has photon-number correlations only with 
itself or with a single matching mode. Such modes, 
which typically are not monochromatic plane waves, have 
been used to describe nonclassical light, mostly for fre¬ 
quency/temporal modes and sometimes also for wavevec- 
tor/spatial modes. Bennink and Boyd introduced the 
term ‘squeezing modes’ for describing squeezing within 
a broad frequency spectrum of a traveling-wave para¬ 


metric amplifier [2]. Opatrny and coworkers used the 
same concept to describe frequency correlations for Kerr- 
squeezed pulses in optical fibres [3]. It is worth stress¬ 
ing that the ‘squeezing modes’ of Ref. [2] as well as the 
‘broadband modes’ of Ref. [3] and the Schmidt modes 
mentioned further here are the same eigenmodes, the 
ones that diagonalize the Hamiltonian producing the ra¬ 
diation. It has been shown [2] that without a proper 
selection of such modes, the degree of measured squeez¬ 
ing decreases considerably. This is a consequence of the 
fragility of squeezing to losses. In the case of multimode 
bright squeezed vacuum (BSV), generated through high- 
gain parametric down-conversion (PDC) [4] or four-wave 
mixing (FWM) [5], the necessity to filter out a single 
Schmidt mode in a clean way, without losing its photons 
or admixing photons from different modes, is especially 
important due to strong thermal fluctuations within each 
mode ms]. If unmatched modes are detected, these fluc¬ 
tuations completely destroy the squeezing. 

There exists a partial solution in the case of homodyne 
detection, where one detects only the modes matching 
the ones of the local oscillator. It is therefore possible 
to select proper modes by tailoring the radiation of the 
local oscillator. In a series of works [9], Fabre, Treps and 
coworkers studied BSV generated by frequency combs 
and its eigenmodes (referred to there as ‘supermodes’). 
They detected the eigenmodes selectively using a spe¬ 
cially tailored local oscillator. Similar strategy of tai¬ 
loring the local oscillator is applied in experiments with 
spatially multimode BSV produced via FWM m- How¬ 
ever, in homodyne detection one cannot make use of the 
radiation in the selected mode, for instance, by coupling 
it with atoms or mechanical systems. If we seek any fur¬ 
ther use of nonclassical light, we need another type of 
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projective filtering, a non-destructive one. 

As a solution, here we consider filtering of a spatial 
mode with an optical fibre. We will show in Section [n| 
that an optical fibre performs the projection of the input 
spatial spectrum on its eigenmode, which can be con¬ 
sidered, to a good approximation, as a Gaussian beam. 
We should stress that only in the case where the fibre 
eigenmode coincides with the radiation Schmidt mode, 
the filtering procedure will retain all features of the ini¬ 
tial radiation such as coherence, peculiarities of the pho¬ 
ton statistics, non-classicality, etc. We test the quality of 
such filtering by measuring the fraction of intensity trans¬ 
mitted through the fibre and comparing it with the the¬ 
oretical expectation, which is the weight of the strongest 
mode in the Schmidt decomposition. As the radiation to 
be tested, we choose multimode BSV generated through 
high-gain PDC. 

It is worth mentioning that the quality of filtering a 
single coherent mode of PDC with an optical fibre has 
been never tested in experiment before. Many authors 
aimed at maximizing the total efficiency of coupling low- 
gain PDC (SPDC) radiation into an optical fibre [11]. It 
has been shown that the optimal case is the one of rather 
tightly focused pump, when SPDC contains few, almost a 
single, mode; but even under this condition the maximum 
coupling efficiency does not exceed 75% m- Besides, as 
we will show further, namely in this case the coupling of 
a single eigenmode is lossy. Other authors [T211T4] max¬ 
imize the heralding efficiency of signal and idler SPDC 
photons coupled into single-mode fibers. It was found 
that the heralding efficiency is maximal and nearly 100% 
for a softly focused pump and, correspondingly, for spa¬ 
tially multimode SPDC m- At the same time, to the 
best of our knowledge, no one considered the losses within 
a single eigenmode. In order to determine the shapes and 
weights of the eigenmodes for PDC radiation, we apply 
the Schmidt decomposition, which is considered in Sec¬ 
tion III We describe the experiment in Section IV where 
we assess the quality of the filtering. 

In the experiment, we deal with spatial (near-held) 
and angular (far-held) modes. Being parts of Fourier- 
conjugate spectra, they are alternative ways to describe 
the radiation. In Section [V| we consider the analogy of 
projective hltering for temporal/frequency spectra and 
propose a method for linear projective hltering of these 
modes. Finally, in Section IVJ we summarize the results. 


II. PROJECTIVE FILTERING AND SPATIAL 
RADIATION EIGENMODES 

A. Filtering with an aperture 

If one needs to select a single mode from the angu¬ 
lar (spatial) spectrum, the simplest strategy is to put an 
aperture of a certain size into the far (near) held. For a 
very small aperture, the mode selected this way will be 
a spherical wave in the hrst case and a plane wave in the 



FIG. 1. Filtering a single eigenmode of broadband PDC ra¬ 
diation with a single-mode hbre. 


second one. It is well known that the radiation after such 
hltering will be coherent (see, for instance. Ref. [T]). Asa 
consequence, the statistical properties (such as photon- 
number distribution) of certain radiation types will be 
also maintained. This will be the case, for instance, for 
thermal light [T]. However, for more fragile types of light, 
e.g. squeezed vacuum, photon-number correlations will 
be lost Moreover, even if the aperture has a spe¬ 

cially chosen size and intensity transmission repeating 
the intensity distribution of a single mode, the hltering 
will still be lossy. 


B. Single-mode fibre as a projective filter 

In order to maintain photon-number correlations, we 
require a diherent strategy, one in which only a single 
eigenmode is hltered out, or a pair of conjugated eigen¬ 
modes, in the case of light with bipartite correlations. 
It is important for the hltering to be of projective type: 
a single held mode of the incident radiation should be 
projected on the eigenmode of the hltering device. This 
kind of hltering is provided by a single-mode hbre m, 
with the restriction that its eigenmode is a Bessel func¬ 
tion, very close to a Gaussian. If a mode of any other 
shape has to be hltered, the hbre could be preceded by a 
spatial light modulator (SLM) performing the transfor¬ 
mation from this shape to a Gaussian [151 US] • The only 
losses introduced this way will be the ones associated 
with the SLM or any alternative device. 

As mentioned before, here we consider PDC as a source 
of multimode radiation (Fig. The PDC radiation 
created in a nonlinear crystal has its near-held ehec- 
tive diameter related to the full width at half maximum 
(FWHM) a of the Gaussian pump. However, as it is 
multimode, its angular divergence is much larger than ex¬ 
pected for a Gaussian beam of waist a. A single-mode h- 
bre hlters out the Gaussian Schmidt mode with the waist 
Wsch losslessly and blocks all other modes. 

Mathematically, the ehect of the hbre is described as 
a projection operation. Let the eigenmode of the hbre in 
the near held be /(r), f |/(r)pdr = 1, with r being the 
coordinate at the input facet of the hbre, and the radi¬ 
ation eigenmodes be J (r)i4/c(r)dr = 6^^- Then 

we can write the photon annihilation operator A in the 
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The modes Un{v) are called coherent modes because the 
radiation within each such a mode is coherent. 

Whenever photon correlations between two beams are 
of interest, especially nonclassical photon-number corre¬ 
lations, rather than coherence within a single beam, Mer¬ 
cer decomposition is not sufficient and the Schmidt de¬ 
composition should be used, as will be shown in the next 
section addressing high-gain PDC. 


FIG. 2. Analogy between multiple orthogonal modes pro¬ 
jected on the eigenmode of a fiber and multiple input modes 
of several beamsplitters projected on a single output mode. 

fibre mode as a linear combination of the photon annihi¬ 
lation operators Ak acting in the radiation eigenmodes, 
which form an ortho normal set: 

A = ^ CkAk. ( 1 ) 

k 

Here, 

Cfe = / f{r)ul{r)dr, y]|C'fc|^ = l, (2) 

k 

are the projections of the fibre mode on the radiation 
eigenmodes. If a single radiation eigenmode rto(r) coin¬ 
cides with the fibre mode, its projection Cq = 1 and the 
other projections are zero. 

Equivalently, the modes can be described in the far 
field as rt/c(q), by introducing the transverse wavevector 
q. Note that the same relation (??) will be valid for 
classical fields instead of the annihilation operators. 

Relation (??) is the same as the one describing the field 
or operator at one output of a series of beamsplitters in 
terms of the fields (operators) at their inputs (Fig. [^: 

n n 

^out ~ ^ ^ ^ ^ \^k\ ~ 1- (3) 

k=0 k=0 

The projections Ck depend then on the field transmis¬ 
sion/reflection coefficients tk^Vk : 

Co = tit2 ...tn,Ci= rit2 ...tn,...,Cn=rn- (4) 


C. Spatial radiation modes 

In the classical case, the eigenmodes of free-space 
radiation are given by the Mercer expansion, defined 
through the first-order Glauber’s correlation function 

min], 

G(i)(ri,r2) = u*(ri)u„(r2). 


III. HIGH-GAIN PDC AND ITS EIGENMODES 

High-gain PDC is a convenient way to produce bright 
squeezed vacuum, a macroscopic quantum state of light 
that is among the most promising sources for quan¬ 
tum technologies. It manifests polarization and photon- 
number entanglement [musi and can violate Bell’s in¬ 
equality under certain experimental conditions [20]. It 
already found applications in quantum imaging m and 
quantum metrology [22], in particular enabling phase 
super-sensitivity [23]. BSV is multimode in angle and 
frequency and can have the mean photon number per 
mode as high as 10 13 [21]. These features provide its 
high information capacity as quantum information can 
be encoded in the photon number of each mode. Ideally 
one would like to isolate and efficiently control each mode 
without losing its nonclassical correlations. 

A. The Schmidt decomposition and the 
Bloch-Messiah reduction 

The eigenmodes of BSV are found from the Schmidt 
decomposition, in which each mode of the signal beam 
is correlated to a single idler mode nansHis]- In signal 
and idler channels taken separately, the modes found this 
way coincide with the coherent modes from the Mercer 
expansion (??) [iZllSo]. 

It has been shown [29] that BSV exhibits the same 
Schmidt modes as two-photon light generated via low- 
gain PDC under the same experimental geometry. The 
modes are found by diagonalizing the Hamiltonian of 
PDC, 

H = ihT j dclsdqiF{cls,c^i)al^^al^.+h.c., (6) 

where F is the coupling parameter scaling as the pump 
amplitude, q^, q^ the signal and idler transverse wavevec¬ 
tor components, the photon creation operators in 

the corresponding plane-wave modes, and F(q 5 ,qi) the 
two-photon amplitude (TPA). This term, as it is used 
here, is conventional: at strong pumping, photons are 
no more emitted in pairs but in large even numbers. The 
Hamiltonian is diagonalized through the Schmidt decom¬ 
position of the TPA, 

= E 'A^^mn{^s)Vmn{q.i), (7) 

m,n 


n 


( 5 ) 
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where are the Schmidt eigenvalues, ^ Xmn = 1, 
and Umni^s)^'^mni^i) ^rc the 2D Schmidt modes of the 
signal and idler radiation. In the degenerate case, where 
signal and idler beams are indistinguishable, r^rnn(q) = 
'^mn (q) • 

After substituting Eq. (??) into Eq. (??), the Hamil¬ 
tonian can be written in the diagonal form (the Bloch- 
Messiah reduction) [29] , 


H = ihV ^ , (8) 

k 


where k = {m,n} and bI are photon creation oper¬ 
ators for the Schmidt modes defined as 

ife = / dqsMfe(q5)aij^, 

f ( 9 ) 

Bi= dqiVkiqi)al., 

In order to label the photons as signal/idler, we assume 
here some additional degree of freedom, for instance, dif¬ 
ferent wavelengths as in non-degenerate PDC. The so¬ 
lution to the Heisenberg equations for the operators in 
Schmidt modes leads to Bogolyubov-type transforma¬ 
tions between the input operators Ako^Bko and the out¬ 
put ones A/e, B^ [29] , 


A/c —C/eA/eO T S]^Bp,Q^ 

Bk =CkBko + 


( 10 ) 


where 


Ck = cosh 




( 11 ) 


and G = jTdt is the parametric gain. 

The mean photon number in mode k is Nj^ = 5^, 
and the total mean photon number in each (signal/idler) 
beam is given by the incoherent sum over separate 
Schmidt modes: N = Thus, the Schmidt eigen¬ 

values at high gain are renormalized [29] . 

_ 4 _ sinh" (Gy^) 


SO that the eigenvalue of a certain Schmidt mode 
determines the contribution of this mode to the total 
number of photons, Nk = N\'^. The effective number 
of Schmidt modes, in the low-gain regime given by the 
Schmidt number K = , is hence reduced at 

high gain, K' = Ej, 

The signal photon annihilation operator A after the fi¬ 
bre is then given by Eq. (??), with Ak being the photon 
annihilation operators in the signal Schmidt modes and 
the projections Ck given by Eq. (??). A similar expres¬ 
sion is valid for the idler photon annihilation operator 
B after the fibre: B = DkBk^ with Dk denoting the 
projections of the idler Schmidt modes on the fibre eigen¬ 
mode. We assume here that the fiber mode function is 


identical for the signal and idler photons. The mean pho¬ 
ton numbers after the fibre are 


TV, = (iti) = ^ \C,f si =Nj2\Ckf K, 

k k 

Ni = (B^B) = ^ \Dkf si = NY,\Dkf X'k- 


(13) 


B. The two-crystal scheme 

In our experiment, described further in Sec. [TV| PDC 
is generated in two consecutive nonlinear crystals placed 
into a common Gaussian pump beam. Eor two crys¬ 
tals of length L at a distance /, the normalized TPA for 
frequency-degenerate or nearly degenerate type-I PDC 
has the form El EH 


P{<ls,q.i) = exp 


X sine 


+ ^i) + cos(^s 

fAk^L\ fAk^L + Ak'l\ 

(—j-A—^ 




X exp{—iAkzL) exp 


{-iAk'l\ 


(14) 

where k^ is the length of signal and idler wavevectors, 
a.re their spherical angles, a is the full width 
at half maximum (EWHM) of the pump intensity dis¬ 
tribution, /Xkz = kp — ko[cos{0s) + cos(6>A] is the lon¬ 
gitudinal mismatch inside each crystal and kp is the 
length of the pump wavevector. In its turn, Ak'^ = 
—/co*’^[cos(0s)+cos(0A] is the longitudinal mismatch 
in the air gap between the crystals, where the pump and 
signal/idler wavevectors take values kp'^'^ and respec¬ 
tively; &si = and no, ng*’^ are the signal/idler 

refractive indices inside the crystals and inside the air 
gap, respectively. 

The Schmidt decomposition is most conveniently found 
in the cylindrical frame of reference, in which the trans¬ 
verse wavevectors are given by their modules, qs,i = 
ko sin Os ,and the azimuthal angles, ,i miEn]. In this 
case, F{q_s^ q^) can be written as a Eourier expansion due 
to its periodicity in (0^ — 0i) [16] . 


n 


where XniQs^Qi) can be found using the inverse Eourier 
transformation. Then, the Schmidt decomposition of 
Xn{qs,qi) yields 


Xn{qsj qi) — ^ ; \/Xmn 

m 


'O^mniQs) 


(16) 


with the functions Umn{Qs) and Vmn{Qi) obeying the nor¬ 
malization condition. 


f 


dqs^mn {Qs)'dkn ) 


f 


dqiVrnniji'i)^ kn^d.i) — ^Ynk' 
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From Eq. (??), we obtain the Schmidt decomposition of 
the TPA (??) with itm„(qs) = = 

Vmn {Qz) „ — in<pi 


with the equality valid only in the case where the Schmidt 
mode and the fiber eigenmode coincide up to a constant 
phase factor. Whenever |Coo| < 1, the filtering procedure 
leads to intrinsic losses. Now let us consider under what 


C. Filtering the first Schmidt mode 

The first Schmidt mode uqq of the BSV state, a Gaus¬ 
sian of waist Wsch^ can be filtered by projecting the angu¬ 
lar spectrum on the eigenmode of a fiber, which is close to 
a Gaussian of waist re, /(q) = exp(—g^/2re^). 

The projections of different Schmidt modes on the eigen¬ 
mode of the fibre are (see Eq. ??) 

p27T pOO 

Cmn= / dtp qdqf{q,4>)umn{q,4>), (17) 

Jo Jo 

Sm n \^mn\‘^ = 1- Notc that siucc the fiber mode does 
not depend on 0, only Schmidt modes with n = 0 will 
contribute in the coupling efficiency. 

Erom (??), the coupling efficiency T for the signal ra¬ 
diation is obtained by summing the photon numbers that 
couple into the fibre from different modes: 

T = (18) 

m,n 

The last inequality follows from the fact that ^>o ^ 
Aqo and shows that the photon number transmitted 
through the fibre is maximized when the fibre mode ex¬ 
actly matches the first Schmidt mode, |(7oo| = 1, and 
Tmax = Aqo- In this case the first Schmidt mode is fil¬ 
tered perfectly, without losses. Only in this case the fil¬ 
tered Schmidt mode maintains all specific features of the 
incoming PDG radiation including non-classicality. 


D. The absence of losses 

However, for the first Schmidt mode to coincide with 
the fibre eigenmode, the former should be real or at least 
have no spatially varying phase, as the latter is a real 
Gaussian function. At the same time, Eq. (??) for the 
TPA includes some phase factors; moreover, even the 
TPA for a single crystal is complex [29] . As a consequence 
the Schmidt modes of down-converted radiation can be 
complex functions. In this case the projection amplitudes 
Cmn entering Eq. (5) of the main text are reduced. Using 
the Gauchy-Schwarz inequality one can show that the 
maximal projection amplitude for the first Schmidt mode 
1^00 0) and the fiber eigenmode /(g, 0) will be still less 

than 1 if the Schmidt mode is complex. Indeed, 

p2tt poo 

ICooP < / d4> qdq\f{q,4>)\'^ X 

Jo Jo 

p2tt poo 

X d4) j qdq\uQo{q,4))\‘^ = 1, (19) 

Jo Jo 


Number of transversal modes 

10 100 



FIG. 3. Numerical estimation of losses accompanying the fil¬ 
tering of PDC generated in the two-crystal system. The pump 
diameter a is varied in order to change the initial number K 
of the transversal modes. The dashed line corresponds to 
the parameters used in our experiment. The insets show the 
intensity distributions of the first Schmidt mode for various 
points in the transverse-wavevector space. 


circumstances the first Schmidt mode would be strongly 
complex. The TPA (??) is complex due to the phase 
factor exj[){—iAkzL — iAk'Jl2}^ which does not depend 
on the pump beam waist a. At the same time, the size of 
the first Schmidt mode does depend on a. By choosing 
appropriate a, one can find the conditions under which 
the phase of the TPA does not vary noticeably within the 
scale of the first Schmidt mode. In this case, the left part 
of (19) will be close to 1 and the losses will be negligibly 
small. 

To demonstrate this, we have performed the Schmidt 
decomposition and calculated Coo foi* various diameters a 
of the pump. For each case, the shape of the first Schmidt 
mode was calculated (the corresponding intensity distri¬ 
butions in transverse wavevector space are shown as in¬ 
sets in Fig. [^. Further, the losses accompanying the 
filtering, given by 1 — |CooP, were numerically calculated 
(shown in Fig. as points connected by a line.) The re¬ 
sult is that the weaker the focusing, the smaller the losses. 
From the inset, one can also notice that for weaker focus¬ 
ing, the first Schmidt mode is narrower and this is why 
the TPA phase in its vicinity is fiat. For focusing into 
more than 100 /rm, the losses are less than 1.5%. In our 
experiment, we use a = 110 /im, which corresponds to 
2% of losses (shown in Fig. [^by dashed line). 

Thus, for a weakly focused pump, the first Schmidt 
mode has a fiat phase and can be filtered nearly loss- 
lessly using a single-mode fibre. In the opposite case of a 
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tightly focused pump, projective filtering with a single¬ 
mode fibre is lossy. In principle, even in this case, the 
phase of the Schmidt mode to be filtered can be elim¬ 
inated before coupling the light into the fibre, but this 
requires some special efforts like using an SLM. 

Note that experiments with SPDC report high (up 
to 96%) heralding efficiencies with a softly focused 
pump despite a low generation rate of photons. 

This is in agreement with our observation (Fig. |^: if 
the heralded mode coincides with the Schmidt mode of a 
highly multimode state, the losses of its fibre coupling can 
be indeed negligible. In contrast, with the pump tightly 
focused into the crystal, so that the resulting number of 
modes is small m, the coupling efficiency is never high, 
again in agreement with our calculations (Fig.[^. At the 
same time, as we have mentioned above, the efficiency of 
filtering a single mode has been never directly measured 
before. The next section considers this measurement. 



IV. EXPERIMENT 
A. Intensity measurements 

The BSV state is generated via high-gain PDC (Fig.[^ 
in two 1 mm long beta barium borate crystals (BBOi, 
BBO 2 ) cut for type-I collinear degenerate phase match¬ 
ing and arranged in the anisotropy compensating con¬ 
figuration at the closest achievable distance of 2.5 ± 0.5 
mm [32113^ . The crystals are pumped by a NdiYAG laser 
third harmonic (wavelength 355 nm, pulse duration 18 ps, 
and repetition rate 1 kHz). The pump is mode-cleaned 
by means of a diamond pinhole (AP) and its polariza¬ 
tion is set by a Gian prism (GP) and a half-wave plate 
HWPi. The pump waist is imaged by lens Lp (100 mm 
focal length) onto the plane between the crystals, where 
its FWHM is 110 um ± 5 um. This experimental config¬ 
uration in the low-gain regime creates a two-photon state 
with the spatial Schmidt number K = 89 (Fig. [^). The 
first Schmidt mode is very close to a Gaussian function, 
with the divergence 0.006 rad and the waist 53 fim. Our 
experiment is performed in the high-gain regime (pump 
power 70 mW), which, according to Eq. (??), leads to the 
reduction in the number of spatial modes. The resulting 
Schmidt number is calculated to be K' = 14.7 (Fig. [^) 
for G = 22.8. 

After the crystals, the pump radiation is cut off by 
a dichroic mirror (DM) and a red glass filter RG 650 
(RG). The spatial filtering is performed by a standard 
step index single-mode fiber (SMF 2 , Thorlabs SM600). 
In the experiment we filter from the PDG radiation, 
with the help of the fibre, Gaussian beams with different 
waists/divergences and measure the fraction T of PDG 
intensity transmitted through the fibre. According to 
Eq. (??), the largest T should equal the first Schmidt 
eigenvalue and should be achieved when the coupled 
Gaussian mode coincides with the first Schmidt mode. 

Eor emulating Gaussian beams with various waists. 


FIG. 4. Setup scheme. The third harmonic of a Nd:YAG 
laser, mode cleaned and properly polarized (aperture AP, 
half-wave plate HWPi and Gian prism GP), focused into crys¬ 
tals BBOi and BBO 2 by lens Lp, generates PDC. An auxiliary 
laser diode beam, mode cleaned by a single-mode fiber SMFi, 
serves to emulate the spatial distribution of different Gaus¬ 
sian modes, with the waists in the same plane. These modes 
are prepared with the help of different lenses Lfc, placed at 
the planes marked by the vertical dotted lines one at a time. 
A Gaussian mode is filtered by the single-mode fiber SMF 2 
after lenses Lc,Lin 2 - Flipping mirrors FMi and FM 2 separate 
the filtering path from the free-space path (dashed lines). Be¬ 
fore PDC enters the monochromator, its polarization is ad¬ 
justed by half-wave and quarter-wave plates (HWP 2 ,QWP 2 ). 
Finally, the signals are analyzed in a Hanbury Brown-Twiss 
interferometer including detectors Di and D 2 . 
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FIG. 5. The Schmidt eigenvalues for PDC from the two- 
crystal scheme used in our experiment (a) in the low-gain 
regime, resulting in the effective mode number A = 89, and 
(b) in the high-gain regime, resulting in K' = 14.7. 


we use an additional GW diode laser with the wave¬ 
length 706.5 nm. The diode laser beam is mode-cleaned 
by a single-mode fiber SMFi, collimated by lens Louti^ 
and then overlapped with the pump beam on the crys¬ 
tals after passing through one of the lenses L/^, with 
the other lenses removed from the beam. Each lens is 
aligned so that it forms the beam waist of a certain di¬ 
ameter Wk (measured by a beam profiler) on the crys- 
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FIG. 6. Spatial filtering of PDC radiation: (a) the coupling 
efficiency versus the angular width /\0k of the Gaussian beam 
coupled into the fiber and (b) versus the longitudinal adjust¬ 
ment dz of the coupling system. The bars at points A and B 
show the absolute uncertainty of the coupling efficiency. 


tals. For each waist diameter, the beam is coupled into 
SMF 2 by means of lens Lc (150 mm focal length) and 
aspheric lens (3.3 mm focal length), with the losses 
12 — 18% including 4% reflection at the uncoated input 
facet. For each lens L/e, efficient coupling of the diode 
laser beam into the fiber indicates that the system fil¬ 
ters out a Gaussian mode with a certain waist and 
the corresponding divergence AOk- After this, the diode 
laser is switched off and the spatial filtering is applied 
to the PDC radiation. Light out-coupled from the fiber 
is sent through a monochromator selecting a bandwidth 
of 0.1 nm around the non-degenerate wavelength 708 
nm. Zero-order quarter-wave plate (QWP 2 ) and half¬ 
wave plate (HWP 2 ) optimize the incoming polarization 
on the monochromator and minimize losses. 

For comparing the data in the presence and in the ab¬ 
sence of filtering, a free-space channel is used where the 
PDC radiation is sent to the monochromator directly, 
avoiding the fiber (Fig. |^. Switching between the free- 
space and spatially filtered channels is done with the flip¬ 
ping mirrors FMl and FM2. The efficiency T of coupling 
into the fiber is measured by dividing the sum signal of 
the detectors Di and D 2 in the presence of filtering by 
the one in its absence. 


Since each Gaussian beam is coupled to SMF 2 with 
85 ± 3% coupling efficiency, the PDC coupling is under¬ 
estimated. We quantify these technical losses and cor¬ 
rect for them at each measured point. In the free-space 
configuration we also measure the total number of spa¬ 
tial modes through the normalized second-order inten¬ 
sity correlation function at the non-degenerate wave¬ 
length of 708 nm, making sure that all angular spectrum 
is collected. Then, the measured correlation function de¬ 
pends on the number of modes K' as = 1 + 1/iF' [34] . 
We obtain g^‘^^ = 1.07 ± 0.01, which indicates K' = 
15 ± 2 spatial modes, in agreement with the calculation 


(Fig.|5 


3 )- 


Fig. 1^ shows the coupling efficiencies for each angu¬ 
lar width AOk of the filtered Gaussian mode as well as 
the theoretical curve, with a good agreement between 
the two. Both curves have their maximum at an angu¬ 
lar width of 6 mrad, which is the angular divergence of 


the first Schmidt mode. Moreover, the maximum cou¬ 
pling efficiency achieved experimentally is 0.201 ± 0.005, 
which matches the first Schmidt eigenvalue Aqq = 0.205 
calculated for our BSV state (Fig.|^. From these values, 
the losses of the filtering procedure do not exceed 2%, as 
expected from the calculation (Fig. [^. This proves that 
filtering of the first Schmidt mode of PDC radiation with 
a single-mode fiber is nearly lossless, up to reflections at 
the fiber facets and imperfect coupling. 

Thus, the efficiency of coupling PDC radiation into the 
fiber is maximal when the first Schmidt mode coincides 
with the fiber eigenmode. One can guess that the first 
Schmidt mode can be targeted by simply maximizing the 
coupling efficiency. In what follows, we show that this 
is indeed the case. In fact, this is the strategy usually 
applied for low-gain PDC m but up to now it has not 
been tested for the filtering of a single mode. 

Starting from a setting where a ‘wrong’ Gaussian mode 
is coupled into the fiber (first point from the right in 
Fig.§i and point A in Fig. [^), we improve the PDC 
coupling by varying the distance between lens Lin 2 and 
the tip of the fiber. This way, we are able to achieve the 
coupling efficiency equal to the first Schmidt eigenvalue 
(point B in Fig. [^). This indicates that the mode col¬ 
lected is indeed the targeted Schmidt mode. Note that 
we achieve this goal by simply moving the fibre tip. At 
first sight impossible, it is feasible in our scheme due 
to the fact that in Gaussian optics, the position of the 
waist image depends on the initial waist size. As a result, 
modification of the diode laser beam waist on the crys¬ 
tals mainly leads to the displacement of the beam waist 
after lens Lin 2 rather than to the change of its size. 


B. Correlation measurements 

One might think that the mode filtering quality could 
be controlled by means of the correlation function mea¬ 
surement [271I2H1I34] since it usually indicates the number 
of modes. However, the measured autocorrelation func¬ 
tion of the signal radiation at zero time delay, 5 fss^( 0 ) 
(further, simply gil\ Fig. [^, is independent of the an¬ 
gular width of the Gaussian mode fully coupled into the 
fiber, and hence of the number of Schmidt modes con¬ 
tributing to the fiber output. 

This can be understood from the analogy with a multi- 
port beamsplitter (Fig.[^. Indeed, with the thermal light 
at the input ports of a beamsplitter, the statistics of light 
at its each output port will be also thermal [35l[36]. Since 
each mode of PDC radiation has thermal statistics, the 
radiation at the output of the fiber will be thermal, and 
one will measure gss^ = 2 regardless of the number of 
modes contributing. Since our frequency filtering is not 

(2) 

perfect, the measured value is lower, gss^ = 1.76T 0.02, 
as shown in Fig. 

This shows that the quality of filtering cannot be as¬ 
sessed from the correlation function measurement for 






FIG. 7. measurement at the non-degenerate wavelength 
^ndeg — 708 um versus the angular width of the Gaussian 
mode coupled into the fiber. 


thermally populated Schmidt modes. At the same time, 
if the modes are populated with two-photon light, corre¬ 
lation measurements will indicate whether one or several 
modes are contributing to the fibre output. 

Indeed, consider the normalized second-order correla¬ 
tion function at zero time delay after the spatial mode 
filter, the auto-correlation function 


(2) {A^A^AA) ^ 

Uss Ar2 

s 

Y, {QnCjrCkCi{AlA]A,Ai), 

^ id,k,l 

and the cross-correlation function 


( 20 ) 


(2) ^ 

“ N,Ni 

^ ^ (QnDjrCkDiiAlBjAkBi). 


( 21 ) 


A straightforward calculation using (??) gives 

— SiCjCf^siSij6}^i -\- SiSjSj^si6i]^6ji. 


Substituting (??) into (??) and (??) we get the following 
expressions for the second-order correlation functions: 


( 2 ) ec:d^s,c. 


(23) 


matches the first Schmidt mode, i.e. \Ck\ = |^/c| = ^ko, 
the expression for simplifies to 

where Nq is the mean photon number in the first Schmidt 
mode. In this case the correlation function after the 
mode filter correctly describes the statistics of a single 
mode of the initial PDC source. However, when the 
mode-matching is not perfect and several modes have sig¬ 
nificantly non-zero coupling coefficients, the correlation 
function value deviates from the expression expected for 
a single-mode PDC field. The dependence on the fiber 



FIG. 8. Dependence of the cross-correlation function g^^^ on 
the size of the mode coupled into the fiber for two values of 
the gain: G = 1 (black line), and G = 10 (red line). The 
width is given in relative units w/wo, where wo is the width 
of the (approximately) Gaussian first Schmidt mode. In both 
cases the Schmidt number A' = 5. 

mode waist, calculated using the double-Gauss model for 
the TPA, is shown in Fig. At first sight, it looks un¬ 
usual that the normalized cross-correlation function is 
minimal for the case of optimal coupling. However, this 
is in line with the usual dependence of normalized cross¬ 
correlation function for PDC on the mean photon num¬ 
ber, given by Eq. (??). At larger photon numbers, PDC 
always has lower normalized cross-correlation function. 
Therefore, optimal coupling, giving the maximal mean 
photon number, at the same time results in a minimum of 
the cross-correlation function. Unfortunately, the depen¬ 
dence flattens out with increasing gain, making the effect 
unobservable under current experimental conditions. 


V. TEMPORAL/FREQUENCY ANALOGY 


The auto-correlation function is identically equal to 2, in 
accordance with the multiport analogy and the results 
of our measurement (Fig[^. The behavior of the cross¬ 
correlation function is more subtle. Its value depends on 
both coupling coefficients between the Schmidt modes 
and the fiber mode and the partial photon numbers in 
each Schmidt mode. In the simplest case where the fiber 


A natural question arises whether a similar technique 
can be applied to filtering a single frequency mode from 
a multimode spectrum. There have been proposals of do¬ 
ing this by means of up-conversion, with the pump mode 
properly tailored in frequency [28l[37]. A ‘quantum phase 
gate’ based on this principle has been recently demon¬ 
strated to provide 80% efficiency for coherent pulses at 












9 


the input [38]. This kind of filtering is also of projective 
type as it projects the field mode on the eigenmode of the 
converter. However, the method is technically compli¬ 
cated, requires phase matching within a broad frequency 
range, may have additional losses for external radiation 
fed into the nonlinear converter and will be also influ¬ 
enced by noise whenever weak input radiation is consid¬ 
ered (a feature of all similar up-converting devices). 

By analogy with the spatial filtering introduced above, 
we can suggest a similar scheme for the filtering of the 
frequency/temporal modes, based on converting the fre¬ 
quency into the angle and then filtering the angle. The 
proposed setup (Fig. is based on a 4f pulse shaper, 
in which the input multimode pulse falls on a diffraction 
grating followed by a cylindrical lens. In the focal plane 
of the lens the vertical coordinate scales linearly with the 
wavelength. In a usual 4f system, a spatially-selective 
device such as a slit or an SLM is placed, after which 
another cylindrical lens and a diffraction grating collect 
the radiation into a single beam. To make the filtering 
projective, one can replace the spatial filter in the middle 
of the 4f scheme by a planar waveguide whose eigenmode 
has a Gaussian profile in the vertical direction. The size 
a of the mode should correspond to the frequency width 
of the Gaussian mode to be filtered. This device will 
project the frequency spectrum on a Gaussian one. If 
a more complicated frequency mode has to be filtered 
out, a mode converter (for instance, an SLM) should be 
placed before the planar waveguide. 



FIG. 9. A setup for linear projective filtering of the frequency 
modes. DG, diffraction grating; CL, cylindrical lens, PW, 
planar waveguide; a, the size of its eigenmode. 


VI. CONCLUSION 

In this paper, we have considered the spatial filtering 
performed by a single-mode fibre and have shown that. 


in contrast to the filtering performed by an aperture, 
it is of projective type and therefore imposes no losses 
on the mode filtered out. An analogy between a fiber 
and a multiport beamsplitter has been drawn. Based on 
this analogy, we have considered this type of filtering for 
the spatial modes of high-gain PDC. It was shown that 
only under the condition that the pump is focused into 
the crystal not too tightly, this filtering can be practi¬ 
cally lossless. Otherwise, if the pump waist is small, the 
Schmidt modes have spatially non-uniform phases and 
the filtering will be lossy unless a phase correction is ap¬ 
plied. 

Further, we have demonstrated projective filtering of 
the first Schmidt mode from the spectrum of high-gain 
PDC in experiment. The total losses accompanying the 
filtering are only caused by imperfect alignment as well as 
the reflection on the uncoated fibre facet and do not ex¬ 
ceed 15%. Importantly, the radiation in the mode filtered 
out this way is not destroyed and is available for further 
use. The method can be extended to higher-order spa¬ 
tial modes by using an appropriate spatial-mode trans¬ 
formations, for instance with the help of a spatial light 
modulator. Furthermore, it has been shown that the cor¬ 
rect Schmidt mode can be filtered simply by maximizing 
the coupling into the fiber, provided that the apertures 
of the lenses do not clip significant portions of the radi¬ 
ation. The transmissivity of the fibre being close to the 
first Schmidt eigenvalue can be therefore used as a cri¬ 
terion for the selection of the first Schmidt mode. On 
the other hand, the autocorrelation function cannot be 
used for this purpose as it is independent on the number 
of modes contributing to the fibre output intensity. The 
cross-correlation function can be used to characterize the 
number of modes contributing but only at low parametric 
gain. 

Finally, a similar technique has been proposed for the 
filtering of a single frequency mode out of the PDG spec¬ 
trum. It is based on a standard 4f pulse-shaping scheme 
where a planar waveguide is used as a spatially selective 
element. 
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